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Abstract. We present a modified description of advection-dominated 
accretion flows with a geometrical viscosity. A simplified treatment of the 
dynamics of the flow in the inner relativistic part is suggested and results 
are compared with the galactic center source Sgr A*. 



1. Introduction 

The existence of a massive black hole in the Galactic Center is a keystone in 
the understanding of stellar dynamics in the inner parsec of the galaxy (see the 
contributions of Ghez and Eckart in this volume) as it is a test case for the the- 
ory of accretion flows in active and non-active galactic nucleii (Mezger, Duschl 
& Zylka, 1996). Standard thin accretion disks fail to explain the luminosity and 
spectrum of the radio source Sgr A* and several models for the accretion and 
emission processes have been suggested. Optically thick synchrotron emission 
from a jet-disk system (Falcke et al., 1996) or from an advection-dominated ac- 
cretion flow (ADAF) were investigated (Narayan et al., 1998), while optically 
thin emission by highly relativistic particles with a quasi-monoenergetic distri- 
bution (e.g. Beckert & Duschl 1997) is also a possible mechanism. Mahadevan 
(1998, and this volume) propose secondary electrons from pion decay in ADAF's 
as the synchrotron emitting particles and showed a way to connect the accretion 
flow as an energy source with the relativistic electrons. We present a simplified 
relativistic treatment of advection-dominated accretion in the galactic center and 
propose a switch between standard and ADAF accretion to distinguish AGN's 
and sub-luminous galactic centers. 

2. Dynamics of Advection-dominated Accretion into Black Holes 

Advection dominated accretion flows are geometrically thick and a description 
in spherical coordinates is appropriate. In order to aviod complications due to 
frame dragging by a rotating black hole we use the Schwarzschild metric and 
hope to see also the characteristics of accretion into a slowly rotating hole. The 
effects of strong gravity and the existence of a horizon are correctly described in 
this way. 
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2.1. Energy Advection and Projectors 



The dynamics of an accretion flow into a black hole is governed by energy and 
momentum conservation. Energy losses due to radiation have to be calculated 
and applied seperately. Angular momentum is transported into the black hole 
and to larger radii than those under consideration so a momentum flux through 
the boudaries of the disk occurs. We will neglect magnetic torques in our treat- 
ment and consider only a contribution of the magentic pressure in the flow to 
the total pressure by scaling it as /3 = PGas/PMag- The dynamics can then be 
derived from the divergence of the stress-energy tensor 



Tg = with T ij = + pg ij - 2pD ij + c- 2 (cpV' + V l Q j ) . (1) 



The density p contains a contribution from the specific enthalpy of the gas like 
p = p{l + with the rest mass density p. V 1 are the components of the 

4-velocity of the flow, p the isotropic pressure and g lJ the metric tensor. The 
viscous part of the stress-energy tensor contains the kinematical viscosity p = up 
and the shear-stresses D ZJ . We have added a contribution by a possible heat 
flux Q l which has to be orthogonal to the velocity. The 0-component of the 
divergence in ([!]) is an energy equation for the flow, but does not resemble the 
second law of thermodynamics. This is formulated for a system at rest so we 
have to project the divergence in ([[]) onto the 4-velocity of the flow. The tool 
for this is a symmetric projector 

= gii + c - 2 v { V j P) (rg) = (2) 

After the projection we get an energy equation like 

pT^- = 2nD ij D i i - Q\ - c~ 2 Q i A i . (3) 

Here we have used the substantial derivative = V^s j and the acceleration 

A % = tj— . In the following we treat stationary and axisymmetric flows with 
radial infall. 

2.2. Isothermal Poloidal Stratification 

In optically thin plasma the heat produced locally by viscous friction will be 
either maintained as internal energy or radiated away. Diffusive transport by 
radiation will not occur. Instead of a radiative transport we expect turbulence 
driven by the magneto-rotational instability (e.g. Balbus & Hawley, 1998) to 
transport heat effectively perpendicular to the flow. So we approximate the 
poloidal structure by an isothermal plasma in hydrostatic equilibrium. The 
poloidal stratification is given by 

|| = p(l + h/c 2 )sm9 cos 9 r 2 Q 2 . (4) 
On 

With the assumed isothermal equation of state we can integrate and obtain 
p(6) = po exp 




(5) 



One obtains the criterion for the flow to be a thin disk as 



Y^T « (^) 2 • (6) 

We take c 2 . = ^ oc T at constant temperature for an ideal gas. The internal 
energy and so the enthalpy rises strongly with temperature when the particles 
become relativistic and the disk thickness will shrink in the relativistic inner 
part close to the black hole. 

2.3. Momentum transport and Viscosity 

In a stationary state as we are interested in the angular momentum distribution 
is determined by diffusive transport due to an effective turbulent viscosity 

dr ^ ^ Sc 2 dr J dr ( dr ) ^ 

For the viscosity we take a parametrisation 

v = a{r-r s )v 4> a ra H' 1 « (8) 

proposed by Duschl, Strittmater & Biermann which is discussed by R. Auer 
(this volume ??). For the self-similiar solution of Narayan &; Yi (1994) this 
viscosity has the same radial dependence as the standard form by Shakura &; 
Sunyaev (1973) . Whenever the ADAF's deviate from the power-law behavior 
the viscosity is not coupled to the thermal state of the flow but rather to the dy- 
namical structure[]. The parameter a is proportional to the inverse of the critical 
Reynolds number 1Z C in the flow. MHD-instabilities may produce turbulence at 
rather small 1Z, so we take a = 0.1. The radial momentum balance 
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rvr)) -[ —vH (9) 



includes the viscous braking in a converging flow and the inertia and pressure 
terms. The centrifugal force reverses sign at r = 3M = 3/211$. The length scale 
M = GA4 /c 2 is the gravitational radius with hA the mass of the black hole. 

2.4. Heating and Cooling of the Flow 

The thermal balance in a stationary state is given by 

Qadv = Qvis — Qheat — Q ■ (10) 



The length scale over which angular momentum is transported by eddie-exchange is limited by 
the distance to the horizon at r s . 
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The flow cools via bremsstrahlung, synchrotron- and inverse Compton losses 
(Q~), while the disk is heated by viscous friction 



m; v 



vYjT 



(ii) 



and the divergence of a possible heat flux Qheat- The heat flux in collisionless 
plasma is approximated by Qheat = ^Scf. In a magnetizied plasma the heat flux 
is expected to be strongly depressed so we take = 10 -3 . The heat generated 
in the flow is not completely radiated but is advected with the flow according to 
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(12) 



The thermodynamic quantities cy , T% are taken from the relativistic equation 
of state for an ideal gas. The thermal balance must be treated for electrons 
and ions separately with an exchange of heat between the two species. This is 
assumed to occur only by Coulomb collisions but is subject to debate, as energy 
transfer due to resonant scattering of electrons by whistler waves can not be 
excluded. 



3. Numerical Treatment 

The vertically averaged dynamical equations for stationary flows presented above 
allow for a self-similiar solution in the Newtonian limit and a constant advection 
ratio (Narayan & Yi, 1994), so that Qadv °c Qtis replaces the energy equation 
(|i~0|). The self-similiarity implies constant logarithmic gradients for all variables 

dQ /dlnft\ Q . . 

; (13) 



dr \ d In r 

as shown here for the angular velocity. The differential equations are decoupled 
and one obtains a set of algebraic equations, which can be easyly solved. We 
use a logarithmic radial grid and solve the algebraic equations with the full en- 
ergy balance at all points. From the local solutions we recalculate the actual 
logarithmic gradients for a next iteration step. In doing so we avoid inner reg- 
ularity conditions at critical (sonic) points. The steep gradients present in the 
solutions at 3 Schwarzschild radii (Fig. |l|) lead to a slow convergence and the 
solutions shown are not completely relaxed to a stationary state. The cooling of 
the flow by radiation is unimportant at very small accretion rates and the radi- 
ation processes have no backreaction on the dynamics of the flow. We include 
bremsstrahlung, synchrotron and inverse Compton-scattered synchrotron radi- 
ation in a rough approximation in the energy equation and derive the spectra 
from the ADAF-solutions afterwards. 



4. Mass Accretion in the Galactic Center 

The accretion rate necessary to power the galactic center source Sgr A* is 
^q-9...-io M Q /y r f or a thi n accretion disk. But upper limits from IR-measure- 
ments allow only a mass flux of 10 -11 M /yr or lower. The inefficiency of 
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Figure 1. Electron (lower curve) and ion temperatures (upper curve) 
as a function of radius and the surface density of an ADAF with a 
mass flow rate of 10~ 5 M Q /yr. Density and temperatures deviate sig- 
nificantly from the self-similiar form at the last stable circular orbit 
at 3Rs where angular momentum does not provide a potential barrier 
anymore. 




Figure 2. Specific Bernoulli number and radial Mach number of the 
flow. The Bernoulli number includes the rest mass of the particles, so a 
zero-energy flow has a Bernoulli number of 1. The Mach number shows 
a transition from subsonic to supersonic behavior inside the last stable 
circular particle orbit. 
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Figure 3. Local luminosity of ADAF and standard thin accretion 
disk for a mass accretion rate of 10~ 5 M /yr. The thin disk luminosity 
(dashed line) is considerably larger than the ADAF emission (solid line) 
and ends at 1.5 black hole radii, where the centrifugal force reverses ist 
direction and is unable to support a stationary thin flow. 



6 



v [GHz] 



Figure 4. A comparision between minimal ADAF model and radio- 
IR observations. We have combined single measurements with variabil- 
ity ranges marked by symbols with circels at both ends. The emission 
by the ADAF itself peaks at 10 3 GHz. Optically thin synchrotron 
emission from relativistic electrons provide the flux at almost all other 
frequencies. The combined spectrum is also shown and falls a little 
short for the measured flux at 100-1000 GHz. 

ADAF-emission processes point to much larger accretion rates, which are deter- 
mined from a simultaneous fit of sub-mm and X-ray measurements. The critical 
parameters are the mass of the black hole, here M = 2.6 10 6 M , the distance 
to the galactic center D = 8 kpc and the hydrogen column density for soft X- 
ray absorption Nn = 6 10 21 cm~ 2 . We assume that most radio observations are 
dominated by optically thin synchrotron emission from relativistic electrons and 
the ADAF contributes significantly in the sub-mm regime only. This minimal 
ADAF-model is calculated with a plasma-/3 = 10 and is correspondingly gas 
pressure dominated. The derived mass accretion rate is 1 • 10 -5 M /yr. Tem- 
peratures and surface density of this model is shown in Fig. |l[ In contradiction 
to the analytical solutions and Narayan & Yi (1995) we find Bernoulli numbers 
less than 1, corresponding to the rest mass at infinity as seen in Fig. |2[ The 
radiation inefficiency is shown in Fig. [|, where the local luminosity is compared 
to the luminosity of a thin disk at the same accretion rate. 



7 



io 8 




Figure 5. The radio to X-ray emission of Sgr A* and the mini- 
mal ADAF-model. ROSAT, Einstein and GINGA observations are 
included. Both the uncorrected and the corrected fluxes corresponding 
to a column density of Ah = 6 10 21 cm~ 2 are shown. The intrinsic 
spectrum of free-free emission is given by the dashed line. According 
to our model most radio observations see optically thin synchrotron 
emission, but the energy output is dominated by the ADAF, both in 
the sub-mm and X-ray regime. The comptonized synchrotron emission 
emerges at IR-optical frequencies and is not included. 
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Figure 6. The total luminosity of standard (upper line) and ADAF's 
as a function of mass accretion rate (left plot) for a 10 8 M© black hole. 
Squares mark calculated luminosities of ADAF's in our treatment. The 
relative luminosity of ADAF's versus M in units of the Eddington 
accretion rate shows the rapid deviation of ADAF's from the standard 
disk luminosity within two decades in the mass flow rate. 

5. Synchrotron and X-ray Spectrum of Sgr A* 

The optically thick self-absorbed synchrotron radiation of the ADAF emerges in 
the radio and sub-mm regime and rises steeply to a maximum at 10 3 GHz. It 
shows an even steeper cut-off for larger frequencies seen in Fig. § Most of the 
radio observations can not be explained by the ADAF itself. One possible source 
for this radiation are relativistic electrons of E = 203 MeV emitting optically 
thin synchrotron radiation in a weak magnetic field B = 2.8 G. Assuming a 
homgeneous source its radius is at least R = 4.5 10 13 cm, much larger than the 
synchrotron emitting region of the ADAF and thus providing an extended halo. 
These relativistic electrons can be accelerated in magnetic reconnection sheets 
(Schopper, Lesch & Birk, 1998) and maintained in a quasi-monoenergetic dis- 
tribution (Jauch, this volume) or they may be secondary electrons from charged 
pion decays (Mahadevan, 1998). 

6. How to switch on an AGN and how to switch it off again 

For mass accretion rates below 0.1 Mead there are two possible states of accretion 
into black holes. The standard thin disk (Frank, King & Raine, 1992 ), which 
looses all the energy produced by viscous friction locally, and the advection- 
dominated accretion flow (ADAF) described above. When we compare the ra- 
dially integrated luminosity of the two kinds of accretion flows at the same 
accretion rate, we find roughly the same luminosity for a critical mass flow rate 
of a 2 MEdd- At smaller rates the luminosity of ADAF's drop dramatically, while 
the luminosity of efficiently radiating flows is proportional to the mass accretion 
rate Lq = 7]Mc 2 . The ratio of the luminosities for ADAF's to Lq gives the 
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amount of specific energy advected into the black hole and the inefficiency of 
the heat transfer between electrons and ions. 

As both states of accretion are available at the same time in a stationary 
flow, there must be a physical reason, which enforces the flow to be advection- 
dominated or not. If a mechanism is at work deciding in which state the flow 
will be, it provides a switch to turn a luminous accreting black hole into a quiet 
sink of matter or the other way round. 

Let us, for the moment being, assume that, whenever available, an accretion 
disk chooses to operate in the advection dominated regime, which is indicated by 
the extremely low luminosity of the Galactic Center. Then comparatively small 
changes in the mass flow rate through the disk may result in dramatic changes 
of the efficiency of the accretion process, i.e., the luminosity of the AGN. This 
highly non-linear behaviour of the luminosity as a function of the accretion rate 
has the potential of acting almost as a switch between an AGN phase of a 
galaxy and a non-active phase whithout having to change the mass flow rate by 
a similar amount as the luminosity. In our example of Fig. ^, at a mass flow 
rate of ~ 1CT 4 M /yr an increase of the rate by two orders of magnitude results 
in almost a million times larger a luminosity. 
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